Single spin detection by qubit SWAP to a molecular nanomagnet 
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Spin state detection is a key but very challenging step for any spin-based solid-state quantum 
computing technology. In fullerene based quantum computer technologies, we here propose to detect 
the single spin inside a fullerene by transferring the quantum information from the endohedral spin 
to the ground states of a molecular nanomagnet Feg, with large spin S=10. We show how to 
perform the required SWAP operation and how to read out the information through state-of-the-art 
techniques such as micro-SQUID. 
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I. INTRODUCTION 

Among the variety of promising technologies to carry 
out quantum information processing, spin-based solid- 
state qubit technologies have great appeal through the 
potential use of large-scale fabrication techniques to scale 
up a quantum computer design Q . Although there have 
been a number of proposals for carrying out universal 
quantum gates in this respect, the read out of the infor- 
mation from individual spin-based qubits - the very last 
and necessary step in any quantum processor routine - is 
still a major difficulty. 

The main idea of current theoretical proposals for sin- 
gle spin detection is to convert spin-state encoded quan- 
tum information into charge state encodings. Since the 
energy of a multi-electron system is spin-dependent, due 
to the Pauli exclusion principle, one can engineer a cur- 
rent whose spin polarization depends on spin-qubit state 
and thusprovides a measurement of the electronic spin 
state 0,U- Experimentally, Scanning- Tunneling- Micro- 
scope electron spin resonance (ESR) has reached the sin- 
gle spin level in the case of iron impurities in Silicon Q . 
However the precise mechanism responsible for this effect 
is not yet well understood. Magnetic Resonance Force 
Microscopy may also become a useful technique for single 
electron spin detection where the cantilever oscillation is 
resonantly driven by a single spin 5] . However this tech- 
nique will require further development before this tech- 
nique reaches single spin sensitivity although significant 
advances have been recently demonstrated [6( . Recently, 
the movement of individual electrons has been evidenced 
by a Cgo-based single electron transistor through electro- 
mechanical coupling 7]. It can be shown that this ar- 
rangement can be used as a single spin detector in the 
presence of an external magnetic field || . 

Although our focus here is to use a Micro-SQUID (Su- 
perconducting Quantum Interference Device - a device 
capable of distinguishing large spin difference (Am = 30) 



we argue that to obtain sufficient sensitivity it will be 
useful to combine this device with methods which could 
convert the small spin-1/2 qubit signal into a system with 
a much larger spin. What we want to detect is the m s 
component of the electronic spin state of the dopant atom 
inside a fullerene Cqq. This doped fullerene or endohe- 
dral fullerene is the primary host for quantum informa- 
tion in a number of fullerene based quantum computer 
designs which include the endohedral molecules N@Ceo 
or P@Ceo jlOL 1 1 it Il2j |. In these types of designs, the 
quantum information is encoded in the nuclear spin of 
the endohedral atom. The execution of quantum gates 
utilises the hyperfine interaction and magnetic dipolar 
coupling via the application of standard nuclear mag- 
netic resonance (NMR) and ESR pulse sequences. It has 
been clearly demonstrated in that fullerene - based 
quantum computing meets the requirement for a quan- 
tum computer, except for the lack of an effective readout 
technique. Since the nuclear spin is less sensitive to the 
external environment, our detection of a single spin state 
would be carried out on the electronic spin state after the 
nuclear spin states have been swapped to the electronic 
spin states. 
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Since the problem we consider is only related to the 
readout of quantum information from a single electron 
spin state, we will neglect nuclear spins in this paper. 
There are three valence electrons in the dopant atom N 
or P, whose ground state in the presence of magnetic 
field is | TIT)- Because other excited states are of much 
higher energy, the only spin state suitable for encoding 
qubits is | ttt), which can be considered as a single spin- 
3/2 state with quartet components \m s — ±3/2) and 
Iro, = ±1/2) in a magnetic field. It has been shown in 
|l2fl that universal quantum information processing can 
be carried out independently with inner qubits | ± 1/2) 
or outer qubits | ± 3/2) . Therefore what we want to 
detect is whether the inner qubit is in |l/2) or | — 1/2) 
or whether the outer qubit is in 1 3/2) or | — 3/2). 



2 



II. COUPLING SYSTEM OF A FULLERENE 

AND A Fes 

Since the electrons of the dopant atoms of C$o are al- 
ways trapped in the cage made of 60 carbon atoms, the 
qubit detection mechanisms suggested in 0, [j| , cannot 
be applied here as these endohedral electrons cannot be 
liberated from the cage without destroying the delicate 
qubit information encoded in their m s spin states. So 
the key step of our scheme is to swap the encapsulated 
spin state to an outside auxiliary spin state which is eas- 
ily detectable. We consider the following Hamiltonian 
describes the endohedral electronic system coupled to a 
nearby auxiliary large-spin system, 



H = H 1 +H 2 + H I 



(1) 



where Hi — —g\[iBB z S\ z denotes the spin (to be de- 
tected), inside the fullerene with gi = 2, \ib being the 
Bohr magneton, B 2 being the magnetic field strength 
along z-axis, and S\ z the Pauli matrix for S — 3/2. H 2 
is related to the auxiliary system with large spin and Hi 
is the coupling between the two systems in a magnetic 
field. 

We employ a molecular nanomagnet made of Fes [13 , 
which is small but has a relatively high magnetic mo- 
ment, as the auxiliary system. It has already been shown 
that one can experimentally prepare a single Fes crys- 
tal in the ground states | ± 10) at very low temperatures 
(< 360 mK) Q. It has also been shown that one can ob- 
serve the quantum tunneling of magnetization in a single 
Fes crystal, which will prove useful for our SWAP oper- 
ation. The lowest energy levels of Fes can be described 
as a spin-10 Hamiltonian 0,0] 



H 2 = -DS% Z 



H tran - g 2 fi B B z S 2z 



(2) 



where D ss 0.275 K is the axial anisotropy constant and 
.92 ~ 5i ■ In what follows, we define u> = gfiBB z in units 
of h = 1, where g = g± ~ g 2 . S 2z is the Pauli matrix 
for S = 10. H tran is the transverse anisotropy due to 
the applied magnetic field in x-y plane, which yields the 
tunnel splitting [lfj . The exact form of H tran is not im- 
portant in our discussion. Moreover, the term concerning 
whose coefficient is much smaller than D is omitted in 
Eq. J2J). We suppose that the detected endohedral spin 
is coupled to the auxiliary spin of the molecular magnet 
by magnetic dipolar interaction, which is generally de- 
scribed as Hi = J (A + B + C + E + F + G) where A = 
(1 - 3 cos 2 9)S lz ®S 2z ,B = -(1/4)(1 - 3 cos 2 9)(S 1+ ® 
S 2 - + Si- <g> S 2+ ) , C = - (3/2) sin cos Oe"** [S u ® S 2+ + 
S 1+ ® S 2z ), E = C*,F = -(3/4)sin 2 0e-* 2 ^ 1+ ® S 2+ , 
and G = -(3/4)sin 2 6»e l2<#, S'i_ ® S 2 - with 6 and the 
usual polar and azimuth angles, respectively. 

Eq. (JIJ can be exactly solved by numerical compu- 
tation. For simplicity, however, we will investigate the 
system under consideration following the assumptions be- 
low: (i) The transverse anisotropy terms will be neglected 
in the following deduction as done in |17j . because it is 



not essential to our conclusion. (We will go beyond this 
assumption when discussing the quantum tunneling of 
magnetization.) (ii) Only the Ising-type term A, related 
to S\ z ® S 2z , will be retained in Hj, which is valid in the 
weak coupling limit 12]. As we will show below, in our 
implementation, the detuning is much larger than Jo, so 
we can safely omit the terms in H/ other than A. As a re- 
sult, the Hamiltonian under consideration can be simply 
rewritten as 

H c = -uj(S 1z + S 2z )- DS\ Z + JS lz <8 S 2z , (3) 

where J = Jo(l — 3 cos 2 9) and we assume J = 0.0175 K 
[l8j . Since all the terms in Eq. © are diagonal, it is easy 
to obtain the eigenstates as well as the corresponding 
eigenenergies, as shown in Table [I] 



III. SWAP IMPLEMENTATION 

Our scheme consists of two essential steps: the first 
step is to swap the detected endohedral electronic spin 
states inside the fullerene to the ground states | ± 10) of 
Fes- Then as the second step, we measure the states 
I ± 10) from Fes- Since a SWAP consists of three 
controlled-NOT (CNOT) operations, we will investigate 
in this section how to carry out CNOT within our model. 
To this end, we will make use of the quantum tunneling 
of magnetization in a single Fes molecule. 

1. Implementation of CNOT 2 \. Tablc^shows the exis- 
tence of degenerate transitions in the subspace regarding 
H 1 . These transitions are heavily dependent on the cou- 
pling (or neighboring) spin states in H 2 . Based on this 
characteristic, the irradiation of a ESR pulse with a reso- 
nant frequency on the first spin yields, in the interaction 
representation, an effective Hamiltonian (h = 1) for the 
subsystem of Hi, 



H — QSix, 
where is the Rabi frequency, 



(4) 
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and the resonant frequency is one of the degenerate tran- 
sition frequencies labeled in Table [I] For a it pulse irra- 
diation of ESR, i.e. fit = it, Hi yields the operator 



P = i 



r 
0010 
0100 

1 0, 



which works independently in the subspace spanned by 
I ± 3/2) or the one spanned by | ± 1/2). 

Therefore, with the operator P, we can flip states 
I ± 3/2) or I ± 1/2) of a single qubit with an ESR pulse 
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whose frequency is determined by the neighboring spin 
state in H 2 - That is actually a CNOT 2 i operation. The- 
oretically, the fidelity of this selective pulse method de- 
pends on more detailed characteristics of the physically 
coupled systems. The implementation time is determined 
by the Rabi frequency f2. 

2. Implementation of CNOTyi- Due to the term S 2z , 
the levels in Fes are not equidistant and so CNOT12 can- 
not be carried out by the above method with selective 
ESR pulses. Let us simply consider the case of very low 
temperature, in which only the ground states | ± 10) 
of H 2 are populated. By considering | ± 10) coupled 
with the possible qubit states of Hi , we show the mag- 
netic field dependence of the energy in Fig. 1, where 
the displayed crossing points correspond to the doublet 
degenerate cases. However, if we include the neglected 
term H tran in our calculation, all of these crossing points 
would turn into avoided crossings due to the tunneling 
splitting. So by sweeping the magnetic field B 2 through 
these (avoided) crossing points, we should have tunnel- 
ing between different magnetization states. For example, 
with the magnetic field B z swept through 0.019 T (from 
0.019_ T to 0.019+ T), we have the magnetization tun- 
neling from |3/2, —10) to |3/2,10). Note that in Fig. 1 
the (avoided) crossing points of different kinds of lines 
mean the places where the magnetization tunneling oc- 
curs with very small probability due to the second or 
higher-order process, e.g. 1 1/2, 10) — * |3/2, — 10) . More- 
over, since the eigenenergies associated with \n, m), n — 
±1/2, ±3/2 and m = ±9, ±8,... are much higher, and 
thus are physically hard to reach in the low temperature 
case, we omit them in Fig. 1. 

The fidelity of our scheme is affected by the follow- 
ing: (1) We must have precise knowledge of the dipolar 
coupling strength J, as this determines the quality of 
the CNOT gates we perform. The J— value can be ob- 
tained experimentally by interrogative ESR pulses. (2) 
We must have very small linewidth selective ESR pulses. 
Since the wavelength of the ESR pulse is much larger 
than the distance between Fes and C60, for avoiding 
exciting Fes when we perform CNOT21, it is required 
that the linewidths of ESR pulses be narrower than J, 
the minimal difference between degenerate transition fre- 
quencies. (3) We must accurately control the tunneling 
time in the performance of CNOT12. This time depends 
on the sweep rate of the magnetic field and the tunnel 
splitting. We prefer a fast tunneling to make the imple- 
mentation time of SWAP shorter than the decoherence 
time. In this case, we need a large tunnel splitting at the 
avoided crossing points. 



IV. DETECTION OF STATES | ± 10) 

Due to the very high sensitivity (which reaches 10 -16 
electromagnetic units |19|L and from the full spectrum 
analysis we have done for nanomagnets |20|. the micro- 
SQUID can hopefully be used to measure the spin states 



of single Fes crystals directly by standard spectroscopy 
with pulsed ESR 0. 



V. DISCUSSION AND CONCLUSION 

Although the experimental value of the coupling J is 
not yet known, we expect J to be suitably large for our 
purpose in implementing the CNOT gates. First, as men- 
tioned above, the prerequisite of a perfect CNOT21 im- 
plementation is that the linewidth of the ESR pulse is 
smaller than J. So the larger the coupling J, the less 
strict the requirement for the linewidth of the ESR pulse. 
Secondly, a larger coupling J is advantageous to achieve 
more accurate implementations of the magnetization tun- 
neling. In terms of our numerical calculation, if J is very 
close to zero, all the (avoided) crossing points would be 
nearly overlapping. It implies that our scheme would 
not work because we could not distinguish different qubit 
spin states from the tunneling signal. On the other hand, 
J cannot be too large. To keep the weak coupling limit in 
Eq. © valid, |J| should be much smaller than min {|w|, 
D}. In our case with J = 0.0175 K corresponding to 350 
MHz, ESR pulses with much narrower linewidths have 
already been achieved experimentally |21|.For the mag- 
netization tunneling |3/2, — 10) <-*■ 1 3/2, 10) occurring at 
B z = 0.019 T, I J| - |w|/2 < D. Nevertheless, as long as 
the dipolar angle 9 in H] is close to tt/2, terms C and E 
in H/ would be nearly zero, and other terms (except A) 
could be neglected in the weak limit and so Eq. © still 
holds. 

An essential point for implementing our scheme is that 
the SWAP time should be shorter than the decoher- 
ence time of the system. For Fes, the Ti times for the 
I ± 10) ground states are very long, while the encap- 
sulated spin states in the fullerene also posses long Ti 
times, 1 sec at 7 K l22jl. a nd probably several seconds 
for lower temperatures |2l|. Thus the dominant source 
of decoherence in our model would be due to the hy- 
perfine level broadening produced by the nuclei [lj, 123 • 
This causes undesired dephasing in our scheme during 
the quantum tunneling of magnetization. We consider 
the SWAP time to be 2tt/VL+ T with T being the mag- 
netization tunneling time. With current experimental 
techniques, f2 = 20 ~ 30 MHz is available [21j, and To 
can be from nanosec to sec depending on the transverse 
magnetic field [24| . In the experiment of [2!| , it is shown 
that, in temperatures below 350 mK, and in the presence 
of B 2 = 42 mT and a transverse B-field 200 mT, the line 
broadening is about 0.8 mT, corresponding to 22.4 MHz 
in units of frequency. So if we assume Q, — 30 MHz, to 
carry out our scheme, we require that To be shorter than 
71 nanosec. 

Another point we should mention is the initial state 
preparation of Fes- Since the detected spin inside the 
fullerene in the readout stage should be well polarized, we 
can simply convert the spin information from H\ to H 2 
by performing CNOT12, instead of SWAP. This would 
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simplify the readout scheme we mentioned above. To 
this end, however, we need to precisely prepare the ini- 
tial states of Fes to be 1 10) or — 110). This too can be 
done by the quantum tunneling of magnetization. When 
we sweep the B z field quickly, we can have an oscillation 
between 1 10} and | — 10), whose frequency is heavily re- 
lated to the tunnel splitting. By stopping the sweep field 
at an exactly chosen time, we can have a perfect initial 
state 1 10) or | - 10) of Fe 8 . 

In summary, we have proposed a potential method to 
efficiently detect the single spin state inside the fullerene 
by means of an auxiliary large spin nanomagnet. Since 
the spin states of the large-spin nanomagnet are measur- 
able with current experimental techniques (e.g. Micro- 



SQUID), we can reliably detect the qubit-encoded spin 
states inside the fullerene by swapping the detected state 
into the ground states of a nanomagnet. 
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TABLE I: The eigenstates and the corresponding eigenener- 
gies obtained from Eq. ©, where the degenerate transition 
frequencies between two nearest neighbor rows in the same 
column are listed in the last row. Due to the axial anisotropy 
term, there is no degenerate transition between two nearest- 
neighbor columns. See text for the details. 
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FIG. 1: Energy versus Bz field plot for the low-lying states at 
very low temperature based on Eq. where D = 0.275A" 
and J = 0.0175-?^. The quantum tunneling of magnetiza- 
tion happens at (avoided) crossing points of the same kinds 
of lines. The (avoided) crossing points of different kinds of 
lines corresponds to the magnetization tunneling with smaller 
probability due to second or higher order process. 



